Gluon distributions in nuclei at small x: guidance from different models 
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The difference between the structure functions measured in nucleons and nuclei 1 1 1 is a very 
important and well known feature of the study of nuclear structure and nuclear collisions. At small 
. . , values of the Bjorken variable x (< 0.01, shadowing region), the structure function F2 per nucleon 

2 \ turns out to be smaller in nuclei than in a free nucleon. This shadowing corresponds to a shadowing of 

Q ■ parton densities in nuclei. While at small x valence quarks are of little importance and the behaviour of 

CN . the sea is expected to follow that of F2, the gluon distribution, which is not an observable quantity, is 

\ badly determined and represents one of the largest uncertainties in computation of cross sections both 

^ ' for moderate and large scales in collineai" factorization |2j. For example, the uncertainty in the 

determination of the glue for Pb at ~ 5 GeV^ at LHC for ?/ = 0, x ~ rriT/y/s ^ 10^^ 10^^, is a 
^ \ factor ~ 3 (see Fig. [0, which for the corresponding cross section in PbPb collisions results in a factor 

~ 9. 

^ I In this situation and while waiting for new experimental data to come from lepton-ion |3l|3|5|| or 

O . pA colliders, the guidance from different theoretical models is of uttermost importance to perform safe 

^ , extrapolations from the region where experimental data exist to those interesting for LHC studies. Two 

different approaches to the problem have been essayed: On the one hand, there exist models which try to 
O ■ explain the origin of shadowing, usually in terms of multiple scattering (in the frame where the nucleus is 

Q . at rest) or parton interactions in the nucleus wave function (in the frame in which the nucleus is moving 



Oh! 



fast). On the other hand, other models parameterize parton densities inside the nucleus at some scale 







p • large enough for perturbative QCD to be appUed reUably, and then evolve these parton densities using 



the DGLAP ||6l0|8l evolution equations; in this way, the origin of the differences of partons densities 
in nucleons with respect to nuclei is not addressed, but contained in the parameterization at which is 
' obtained from a fit to experimental data. 



^ ■ 1. Multiple scattering and saturation models 

;h ■ 

I The nature of shadowing is well understood qualitatively: In the rest frame of the nucleus, the incoming 

photon splits, at small enough x, into a pair long before reaching the nucleus, with a coherence length 
/c oc 1/ {mj\[x) with mj\[ the nucleon mass, which at small enough x becomes of the order of or greater 
than the nuclear size. Thus this qq pair interacts coherently with the nucleus with typical hadronic cross 
sections, which results in absorption ||9l[l0l[ni[l2l[l3l (see llT4l for a simple geometrical approach in this 
framework). In this way nuclear shadowing is a consequence of multiple scattering and is thus related 
with diffraction (see e.g. 031(161 '). 

Multiple scattering is usually formulated in the dipole model fT3llT7l . which is equivalent to kx- 
factorization [18| at leading order. In this framework the 7* -nucleus cross section is expressed through 
the convolution of the probability of the transversal or longitudinal 7* to split into a qq pair of transverse 
dimension r times the cross section for scattering of such dipole with the nucleus. It is this dipole- 
nucleus cross section at fixed impact parameter which saturates (i.e. gets a maximum value allowed by 
unitarity), most frequently by multiple scattering in the Glauber-Gribov approach [19 20 21 16|. This 
dipole-target cross section is related through a Bessel-Fourier transform to the so-called unintegrated 
gluon distribution ipAix, kx) in /ct -factorization l22l[T9ll20l . which in turn can be related to the usual 
collinear gluon density through 

xGa{x,Q^)= dk^ ipA{x,kT) (1) 
Ja? 



(with some infrared cut-off, if required) although this identification is only true for large » 
Q\ 1231 124II . with the saturation momentum corresponding to a transverse length scale where the 
saturation of the dipole-target cross section occurs. 

Other formulations of multiple scattering do not use the dipole formulation but relate shadowing 
with diffraction | iT TF] by Gribov theory. In this way diffraction in lepton scattering on nucleons is 
related to nuclear shadowing. 

An equivalent explanation to multiple scattering in the frame in which the nucleus is moving 
fast, is that gluon recombination due to the overlap of the gluon clouds from different nucleons, makes 
gluon density in nucleus with mass number A smaller than A times that in a free nucleon l25l l26l : 
at small x the interaction develops over longitudinal distances z ~ Xjimf^x) which become of the 
order of or larger than the nuclear size, leading to the overlap of gluon clouds from different nucleons 
located in a transverse area ~ l/Q^- These studies have received great theoretical impulse with the 
development of semiclassical ideas in QCD and the appearance of non-linear equations for evolution 
in X in this framework (see f27','28l and references therein), although saturation appears to be different 
from shadowing |23, 24] (i.e. the reduction in the number of gluons as defined in DIS is not apparent). 
In the semiclassical framework the gluon field at saturation reaches a maximum value and becomes 
proportional to the inverse QCD coupling constant, gluon correlations are absent and a form for the 
dipole-target cross section appears which leads to a geometrical scaling of 7* -target cross section; such 
scaling has been found in small x nucleon data [291, and also in analytical and numerical solutions of 
the non-linear equations l30l 131113211331 . but apparently the region where this scaling should be seen 
in lepton-nucleus collisions has not been reached yet in available experimental data l34ll . Apart from 
eA colliders |3','4','5'|, LHC will be the place to look for non-linear effects. However, the consequences 
of such high density configurations may be masked in AB collisions by other effects, as those due to 
final state interactions. pA collisions would be, thus, essential for this type of studies (see e.g. |35|). 
Moreover, they would be required to fix the baseline for other studies in AB collisions as QGP search 
and characterisation. As a final comment, let us indicate that the non-linear terms in fSS' "Sf) | are of a 
higher-twist nature and in this case the low density limit recovers the DGLAP equations |6 7, 8|, while 
the non-linear equations for evolution in x I27II28II do not correspond to any definite twist and their linear 
limit correspond to the BFKL equation I36ll37l . 

Let us comment a little more on the difference between shadowing for gluons and saturation. In 
saturation models |27 28 1 saturation, defined by a maximum value of the gluon field or by the scattering 
becoming black, and shadowing for gluons, defined by the ratio xGa/ {xGn) being < 1, are apparently 
different phenomena, i.e. saturation does not necessarily lead to shadowing for gluons l23ll24l . Indeed, 
in the framework of numerical studies of the non-linear equations for small x evolution, the unintegrated 
gluon distribution turns out to be a universal function of just one variable t = k'^/Ql I30ll31l l32 1 (in 1 33l 
this universality is analytically shown to be fulfilled up to kj, much larger than Ql), vanishing quickly 
for fcj. > Q^; this scaling appears also in the analysis of DIS experimental data on nucleons l29ll and has 
been searched for in nuclear data |34|. In nuclei, Ql increases with increasing nuclear size, centrality 
and energy. So, through the relation with the collinear gluon given by Eq. Q, this scaling implies that 
the integral gives the same value (up to logarithmic corrections if a perturbative tail oc l/fcy exists for 
large k^^ ^> Ql), provided that the upper Umit of the integration domain Q^ ^ Q^. Thus the ratio 
xGa/{xGn) keeps equal to 1 (or approaching 1 as a ratio of logs if the perturbative tail exits) with x 
for Q^ » Q^. So saturation or non-linear evolution do not automatically guarantee the existence of 
shadowing for gluons. As shadowing is in many models described in terms of multiple scattering, this 
result I23II24II looks somewhat surprising. 

2. DGLAP evolution models 

On the other hand, a different approach is taken in ^38", "39", '40', '4T| : parton densities inside the nucleus 
are parameterized at some scale ~ 1 5 GeV^ and then evolved using the DGLAP l6||7j|8J evolution 



equations. In this way, all nuclear effects on parton densities are included in the parameterization at 
Qq, which is obtained from a fit to experimental data. The differences in the approaches come mainly 
from the sets of data used (e.g. the use of Drell-Yan data or not) to constrain the parton distributions. In 
these calculations the lack of experimental data makes the gluon to be badly constrained at very small 
x; experimental data on the evolution of F2 with logQ^ give direct constrains |42| through DGLAP 
evolution but only for x > 0.01. Please refer to l42ll and references therein for more discussions on this 
kind of models. 

3. Comparison between different approaches 

The results from different models usually depend on additional semiphenomenological assumptions and 
often contradict each other. For example, concerning the -dependence of the effect, in 19., .10. Ill 
(121 [TSll it is argued that qq configurations of a large dimension give the dominant contribution to the 
absorption, which results essentially independent of Q^. On the other hand, in the gluon recombination 
approach of f25"26 1 the absorption is obtained as a clear higher-twist effect dying out at large Q^. Finally, 
in the models |38,^9.,40, 41| which use DGLAP, all -dependence comes from QCD evolution and 
is thus of a logarithmic, leading-twist nature. Predictions (particularly for the gluon density) on the 
x-evolution towards small x turn out to be very different. 

Let us compare between different approaches whose numerical results are available (when con- 
sidering an approach based on the dipole model as in (19], one should keep in mind the difficulties to 
identify at small and moderate the integral of the unintegrated gluon distribution with the ordinary 
gluon density, see Eq. (Q and comments below it). A comparison at = 5 GeV^ for the ratio of 
gluon densities in Pb over proton, can be found in Fig. ^ There it can be seen that the results of 
fT9"38"39"40"/2r| at j; ~ 10^^ (relevant for RHIC) roughly coincide, while they are higher than those 
of L16. 43J ; at x ~ 10^* 10^^ (accessible at LHC) the results of 1 19| become smaller than those of 
l38ll39ll4()ll2T]l . get close to those of (T6i and are still larger than those of |43 1. Apart from the constraints 
coming from existing DIS experimental data on nuclei which are very loose for the glue at small x, in 
f38"39"40'| the saturation of gluon shadowing comes mainly from the initial condition for DGLAP evo- 
lution where this saturation has been imposed, while [ 19. J6. .21 1 are multiple scattering models. In l,43l 
the behaviour of the glue, in the form of a Q^-independent parameterization, has been fixed at x ~ 10~^ 
in order to reproduce charged particle multiplicities in AuAu colUsions at RHIC. In this approach the 
strongest gluon shadowing is obtained. However, it seems to be in disagreement with existing DIS data 
|42|. Additional caution has to be taken to compare results from multiple scattering models with those 
coming from DGLAP analysis l38l l39l l40l : the ratios for the glue at some moderate, fixed and very 
small X may result smaller (e.g. in fT9|) than the ratios for F2 at the same x, Q^, which might lead to 
problems with leading-twist DGLAP evolution, see |42I . 
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Fig. 1: Ratios of gluon distribution functions from different models at Q^=5 GeV^; HKM refers to the results from 1401 . 
Sarcevic to those from 1211 . EKS98 to those from I38II39I . Frankfurt to those from 1 16 1, Armesto to those from |19||20!| and 
new HIJING to those from 1431 . The bands represent the ranges of x relevant for processes with a scale = 5 GeV^ at RHIC 
and LHC. 
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